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^ , Abstract 

' In this paper, we study the orbital stability for a four-parameter family of periodic stationary 

. traveling wave solutions to the generalized Korteweg-de Vries (gKdV) equation 

ON- u t = u xxx + f(u) x . 



PL, 



X 



In particular, we derive sufficient conditions for such a solution to be orbitally stable in terms of 
the Hessian of the classical action of the corresponding traveling wave ordinary differential equation 



■ restricted to the manifold of periodic traveling wave solution. We show this condition is equivalent to 

the solution being spectrally stable with respect to perturbations of the same period in the case when 
/(it) = u 2 (the Korteweg-de Vries equation) and in neighborhoods of the homoclinic and equilibrium 
solutions if /(«) = u p+1 for some p > 1. 



1 Introduction 

> 

This paper concerns the stability analysis of periodic traveling wave solutions of the generalized Korteweg- 
de Vries (gKdV) equation 

^ ' ut = u xxx + f(u) x (1) 

where / is a sufficiently smooth non-linearity satisfying certain convexity assumptions. Probably the 
most famous equation among this family is given by f(u) = u 2 , in which case {1} corresponds to the 
Q^ \ Korteweg-de Vries (KdV) equation. The KdV serves as an approximate description of small amplitude 

waves propagating in a weakly dispersive media. When the amplitude of the wave is not small, however, 
the gKdV (H|) arises where the nonlinearity / describes the effects of the solution interacting with itself 
through the nonlinear media. 

It is well known that the gKdV equation admits traveling wave solutions of the form 



u(x,t) =u c (x + ct), feM (2) 

for wave speeds c > 0. Historically, there has been much interest in the stability of traveling solitary waves 
of the form ^ where the profile u c decays exponentially to zero as its argument becomes unbounded. 
Such waves were initially discovered by Scott Russell in the case of the KdV where the traveling wave is 
termed a soliton. While ([TJ) does not in general posess exact "soliton" solutions, which requires complete 
integrability of the partial differential equation, exponentially decaying traveling wave solutions still exist. 
Moreover, the stability of such solitary waves is well understood and dates back to the pioneering work of 
Benjamin [3], which was then further developed by Bona [5], Grillakis [12], Grillakis, Shatah and Strauss 
[13] , Bona, Souganides and Strauss [B], Pego and Weinstein [TIB US] j Weinstein |221 123j . and many others. 
In this theory, it is shown that traveling solitary waves of ([1]) are orbitally stable if the the solitary wave 
stability index 

d f°° 

- u c dx (3) 
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is positive, and is orbitally unstable if this index is negative. In the case where ([T]) has a power-law 
nonlinearity f(u) = u p+1 , the sign of this stability index is positive if p < 4 and is negative if p > 4. 
Moreover, in [THl HO] it was shown that the mechanism for this instability is as follows: Linearizing the 
traveling wave partial differential equation 

u t = u xxx + f(u) x - cu x , (4) 

which is satisfied by the traveling solitary wave profile, about the solution u c and taking the Fourier 
transform in time leads to a spectral problem of the form 

d x C[u c ]v = fiv 

considered on the real Hilbert space L 2 (M.), where C[u c ] is a second order self adjoint differential operator 
with asymptotically constant coefficients. The authors then make a detailed study of the Evans function 
D(fi), which is an analytic function such that if if -0 is a solution of ((3]) satisfying ip{x) ~ e ux as x — > oo, 
then ip(x) ~ D(fi)e UJX as x — » — oo: in essence, D(n) plays the role of a transmission coefficient familiar 
from quantum scattering theory. This approach motivated by the fact that for Re(/i) > the vanishing 
of D([i) implies that \i is an L 2 eigenvalue of the linearized operator d x C[u c ], and conversely. Pego and 
Weinstein were able to use this machinery to prove that the Evans function satisfies 

lim sign(D(/j)) > 

as well as the asymptotic relation 

DQt) =C 1 (J^ J u c (x) 2 dx^j + o(| M | 2 ) 

in a neighborhood of y, = 0, for some positive constant C\. Thus, in the case when the solitary wave 
stability index is negative, it follows by the continuity of D(fi) for fj, £ R + that D(fi) < for small positive 
\x and hence D(p) must have a positive root, thus proving exponential instability of the underlying 
traveling solitary wave in this case. 

In this paper, however, we are concerned with traveling wave solutions of ([1]) of the form lp]). where 
this time we require the profile u c be a periodic function of its argument. In contrast to the traveling 
solitary wave theory, relatively little is known concerning the stability of periodic traveling waves of 
nonlinear dispersive equations such as the gKdV. Existing results usually come two types: spectral 
stability with respect to localized or bounded perturbations, and orbital (nonlinear) stability with respect 
to periodic perturbations. Most spectral stability results seem to rely on a Floquet-Bloch decomposition 
of the linearized operator and a detailed analysis of the resulting family of spectral problems, or else 
perturbation techniques which analyze modulational instability (spectrum near the origin) . Many of the 
stability results known to the author for the traveling wave solutions of the gKdV concern the stability 
of the cnoidal solutions of the KdV 

u(x, t) = uo + 12k 2 n 2 cn 2 (k (x - x + (8k 2 n 2 - 4k 2 + u ) t) ,k) , 

where k £ [0, 1) and k, xq, and uq are real constants. Such cnoidal solutions of the KdV have been 
studied by McKean [18j , and more recently in papers by Pava, Bona, and Scialom [2] and by Bottman and 
Deconinck [7j . The results in [18] uses the complete integrability of the KdV in his study of the periodic 
initial value problem in order to show nonlinear stability of the cnoidal solutions to perturbations of the 
same period. Also using the machinery of complete integrability, in [7] the spectrum of the linearized 
operator on the Hilbert space L 2 (M.) is explicitly computed and shown to be confined to the imaginary 
axis. In particular, it follows that cnoidal solutions of the KdV are spectrally stable to perturbations of 
the same period, and more generally, perturbations with periods which are integer multiples of the period 
of the cnoidal wave. 
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In [2], more classical energy functional methods are used to show that such cnoidal solutions are 
orbitally stable to perturbations of the same period as the cnoidal wave. Since their techniques do not 
rely on the complete integrability of the KdV, it is possible that their methods can be extended to 
equations in the family described by (|TJ) other than the KdV. In this paper, conduct exactly such an 
analysis and show how the orbital stability of the cnoidal solutions of the KdV can be recovered from 
this theory once it is known that such solutions are spectrally stable to perturbations of the same period. 
It should be noted, however, that our analysis varies slightly from that of [2J due to the fact that the 
§Q is not in general Galilean invariant. Thus, we must make use of all three conserved quantities of the 
gKdV flow: see Remark [7] below. Our goal is to derive geometric criterion for the orbital stability of such 
solutions which are universal to the family of equations ([I}. As such, our results are not as explicit as in 
the papers mentioned previously, yet they apply to a much broader range of equations. 

Returning to the generalized KdV equation |T]) , spectral stability results have recently been obtained 
by Haragus and Kapitula [T^] as well as by the author in collaboration with Bronski [5]. In [TB], the 
spectral stability of small amplitude periodic traveling wave solutions of ([1]) with f(u) = u p+l was 
studied. By using a Floquet-Bloch decomposition of the linearized spectral problem, the authors found 
that such solution^ are spectrally stable if p S [1,2) and exhibit a modulational instability if p > 2. 
In particular, they found that such solutions are always spectrally stable to perturbations of the same 
period: in section 5, we will verify and extend this result through the use of the periodic Evans function. 
In [S] , a modulational instability analysis of the periodic traveling wave solutions of ([I]) by using Floquet 
theory and developing a perturbation theory appropriate to the Jordan structure of the period map at the 
origin. As a by product of their analysis, a sufficient condition for exponential instability of the underlying 
periodic traveling wave with respect to periodic perturbations was obtained in terms of the conserved 
quantities of the gKdV flow. In particular, a stability index was derived in a manner quite similar to the 
solitary wave theory outlined above such that the negativity of this index implies exponential instability 
of the periodic traveling wave with respect to perturbations of the same period. The relevant results of 
this analysis can be found in section 3. 

It seems natural to wonder role this periodic instability index derived in [5] plays in the nonlinear 
stability of the periodic traveling wave. As mentioned above, the analogue of this index controls the 
nonlinear stability in the solitary wave context. Thus, one would like the periodic traveling wave to be 
nonlinearly stable whenever the aforementioned periodic stability index positive. While we are able to 
show this is true in certain cases, we find that two other quantities, which are essentially not present in 
the spectral stability theorjH nor the solitary wave theory, play a role in the nonlinear stability. This is 
the content of our main theorem, which is stated below [j. 

Theorem 1. Let u{x + ct) be a periodic traveling wave solution of ([I]), whose orbit is bounded in phase 
space by the homoclinic orbit. Moreover, assume the periodic instability index {T, M, P} a ,E,c is positive. 
Then there exist two real real numbers, denoted Te and {T, M} at E, associated to the periodic traveling 
wave profile u such that if Te > and {T, M} a ^E > 0, then there exists Cq > and e > such that for 
all 4>q S X with || 0o II -X" < £, the solution 4>{x,t) of {1} with initial data <j)(x,0) — u(x) + (f)(x) satisfies 



Remark 1. Throughout this paper, "orbital stability" will always mean orbital stability with respect to 
perturbations of the same period as the underlying wave. 

1 In the next section we will give a parametrization of the periodic traveling wave solutions of JTJ. In 1161 . they only 
consider small amplitude periodic traveling waves with the integration parameter a (see section 2 for definition) being small. 

2 This is not quite correct. They are present, but their signs do not play into the spectral stability theory. See section 3 
for more details. 

3 For the definition of the real Hilbert space X, see section 4. 




u(x + ct + £)\\x < C o ||0o||x 



for all t > 0. 
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The outline for this paper is as follows. Section 2 will be devoted to a study of the basic properties of 
the periodic traveling wave solution of equation ([I]). In section 3, we will recall the recent results of [8] 
concerning the spectral stability of periodic traveling wave solutions of {1} with respect to perturbations 
of the same period. The resulting instability index will play an important role throughout the rest of 
the paper. Section 4 is devoted to the proof of Theorem [TJ Finally, two applications of our theory 
are described in sections 5 and 6 in the case of a power-law nonlinearity f{u) — u p+l for p > 1. In 
section 5, we study the orbital stability of periodic traveling wave solutions of (JTJ) in neighborhoods of 
the homoclinic and equilibrium solutions. Section 6 is devoted to the application of our theory to the 
case of the KdV. In particular, it is shown that such solutions are orbitally stable if and only if they are 
spectrally stable to perturbations of the same period as the underlying wave. 

2 Properties of the Stationary Periodic Traveling Waves 

In this section, we recall the basic properties of the periodic traveling wave solutions of (p}. For each 
number c > 0, a stationary traveling wave solution of ([I]) with wave speed c is a solution of the traveling 
wave ordinary differential equation 

Uxxx + f(u) x - cu x = 0, (5) 

i.e. they are solutions of ([TJ which are stationary in the moving coordinate frame defined by x + ct. 
Clearly, such solutions are reducible to quadrature and satisfy 

u xx + /(«)- cu = a, (6) 
X -u\ + F(u) - ~u 2 - an = E, (7) 

where a and E are real constants of integration, and F satisfies F' = /, F(Q) = 0. In order to ensure the 
existence of periodic orbits of (JSJ), we require that the effective potential 

V(u; a, c) = F(u) — —u 2 — au 

is of class C 2 (R) and has a local minimum. Notice this places a restriction on the allowable parameter 
regime for our problem. This motivates the following definition. 

Definition 1. We define the set f2 C K 3 to be the open set consisting of all triples (a,E,c) such that 
c > and the solution u(x) = u(x;a,E,c) of (J7|) is periodic and its orbit in phase space lies below the 
separatrix. 

Remark 2. Taking into account the translation invariance of ([T]). it follows that for each (a,E,c) G fl 
we can construct a one-parameter family of periodic traveling wave solutions of |T]) : namely 

u^(x, t) — u(x + ct + £;a,E,c) 

where £ £ R. Thus, the periodic traveling waves of ([T]) constitute a four dimensional manifold of solutions. 
However, outside of the null- direction of the linearized operator which this generates, the added constant 
of integration does not play an important role in our theory. In particular, we can mod out the continuous 
symmetry of ([T]) by reguiring all periodic traveling wave solutions to satisfy the conditions u x (0) = and 
V '(it(0)) < 0. As a result, each periodic solution of is an even function of the variable x with a local 
maximum at x = 0. 

Remark 3. Notice that a and E are conserved guantities of the ODE flow generated by ([5]). Moreover, 
the classical solitary waves corresponding to solutions of ([T]) satisfying lim^^j-oo u(x) = correspond to 
a = E = 0, and hence constitute a codimension two subset of the traveling wave solutions of ([TJ. It 
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seems natural then to expect that the stability of periodic traveling wave solutions will involve variations 
in these added parameters, just as the solitary wave stability index involves variations in the one (modulo 
translation) free parameter c. 

Throughout this paper, we will always assume that our periodic traveling waves correspond to an 
(a, E, c) within the open region £7, and that the roots u± of E = V(u; a, c) with V(u; a, c) < E for 
u G (tt_,U-|_) are simple. Moreover, we assume the potential V does not have a local maximum in the 
open interval u + ). It follows that u± are C 1 functions of a, E, c on fi, and that it(0) = U-. Moreover, 
given (a, E, c) € f2, we define the period of the corresponding solution to be 

r + du 

T = T(a,E,c):=2 . 8 

The above interval can be regularized at the square root branch points u_, u + by the following pro- 
cedural: Write E — V(u;a,c) = (u — u_)(u + — u)Q(u) and consider the change of variables u = 
"++"- 4- "+~"~ sin(fl). Notice that Q(u) ^ on [u_,it + ]. It follows that du — y/(u — — u)dO 

and hence ([5]) can be written in a regularized form as 



T{a,E,c) = 2 J 



d0 



Q «±+^ + H±^ sin ^) 



In particular, we can differentiate the above relation with respect to the parameters a, E, and c within 
the parameter regime Q. Similarly the mass, momentum, and Hamiltonian of the traveling wave are 
given by 



T ru + u du 



o 

,2\ _ / „.2/ 



M(a, E, c) = {u)= I u(x) dx = 2 — — (9) 

J«_ a/2 (E — V(u; a, cj) 



P(a, E, c) = (u J ) = / u\x) dx = 2 = (10) 

Jo Ju- \J2 [E — V(u; a, c)) 

H(a, E,c) = (4- F(u)\ = 2 f + ^zXtMzM du. (11) 
V ' \ 2 V 7 7„_ v/2(£;-V( W ;a,c)) V ; 

Notice that these integrals can be regularized as above, and represent conserved quantities of the gKdV 
flow restricted to the manifold of periodic traveling wave solutions. In particular one can differentiate 
the above expressions with respect to the parameters (a,E,c). 

Throughout this paper, a large role will be played by the gradients of the above conserved quantities. 
However, by the Hamiltonian structure of ([5]), the corresponding derivatives of the period, mass, and 
momentum restricted to a periodic traveling wave u(-;a,E,c) with (a,E,c) € f2 satisfy several useful 
identities. In particular if we define the classical action 

/r T /•«+ 
u x du = u\dx = 2 \ \J2(E - V(u; a, c)) du 
JO Ju- 

(which is not itself conserved) then this quantity satisfies the following relations 

K E (a,E,c)=T{a,E,c) 
K a (a,E,c)=M(a,E,c) 
2K c (a,E,c) = P(a,E,c). 



4 In section 5 we will use a slightly different regularization which is well-adapted to the KdV equation. 
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Moreover, using the fact that T, M, P and H are C 1 functions of parameters (a, E, c), the above implies 
the following relationship between the gradients of the conserved quantities of the gKdV flow restricted 
to the periodic traveling waves: 

£JVT + aVM + -VP + V-ff = 0, (12) 

where V = (<9 a , c>e, <9 C ): see the appendix of [5] for details of this calculation. The subsequent theory is 
developed most naturally in terms of the quantities T, M, and P. However, it is possible to restate our 
results in terms of M, P, and H using the identity (fT2|) . This is desirable since these have a natural 
interpretation as conserved quantities of the partial differential equation (P) . 

We now discuss the parametrization of the periodic solutions of ([5]) in more detail. A major technical 
necessity throughout this paper is that the constants of motion for the PDE flow defined by {1} provide (at 
least locally) a good parametrization for the periodic traveling wave solutions. In particular, we assume for 
a given (a, E, c) € fi the conserved quantities (H, M, P) are good local coordinates for the periodic travel- 
ing waves near (a, E, c). More precisely, we assume the map (a, E, c) — * (H(a, E, c), M(a, E, c),P(a, E, c)) 
has a unique C 1 inverse in a neighborhood of (a, E, c). If we adopt the notation 



{/,<?}* 



fx 9x 

fy 9y 



for 2x2 Jacobians, and {/,<?, h} x . yyZ for the analogous 3x3 Jacobian, it follows this is possible exactly 
when {H, M, P} a ,E,c is non-zero, which is equivalent to {T, M, P} a ,E,c ^ by (fT2"|) . Also, we will need 
to know that two of quantities T, M and P provide a local parametrization for the traveling waves with 
fixed wave speed. By reasoning as above, this happens exactly when the matrix 

T a M a P a 

T E M E P e 

has full rank. A sufficient requirement is thus {T, M} a ^E ^ 0: trivial modifications of our theory is needed 
if {T, M} at E — but one of the other sub-determinants do not vanish. 



3 Spectral Stability Analysis 

In this section, we recall the relevant results of [5] on the spectral stability of periodic-traveling wave 
solutions of the gKdV. Suppose that u = u( ■ ; a, E, c) £ C 3 (IR; R) is a T-periodic periodic solution of ((H). 
Linearizing about this solution and taking the Fourier transform in time leads to the spectral problem 

d x L[u]v = fiv (13) 

considered on L 2 (M; K), where C[u] := —d\ — f'( u ) + c is a closed symmetric linear operator with periodic 
coefficients. In particular, since u is bounded it follows that C[u] is in fact a self-adjoint operator on L 2 (M.) 
with densely defined domain C°°(1R). Notice that considering (|T5|) on L 2 (M.) corresponds to considering 
the spectral stability of u with respect to localized perturbations. One could also study the spectral 
stability with respect to uniformly bounded perturbations, but by standard results in Floquet theory the 
resulting theories are equivalent. In order to understand the nature of the spectrum of (H3J) we make the 
following definition. 

Definition 2. The monodromy matrix M(/i) is defined to be the period map 
where <i>(a;,/z) satisfies the first order system 



= H(rr, M )$, $(0,/i) = 1 



(14) 
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where I is the 3x3 identity matrix and 

/ 10 

H(x,fjt) = 10 

V -/* -««/"(«) -/'(«) +c o 

It follows from a straightforward calculation that the operator has no L 2 eigenvalues. Indeed, 

from Floquet's theorem we know that every solution of (| 14[) can be written in the form <&(a;) = e u)x p(x) 
where p is a T-periodic function and A = e wT is an eigenvalue of M(/i). Conversely, each eigenvalue of 
M(/i) leads to a solution of this form. It then follows that for any N £ N one has 

$(JVT) = e" NT P (0) = X N p(0). 

If $(x) decays as a; — ► oo, the above relation implies that it must be unbounded as x — > — oo, from which 
our claim follows. In particular, a solution $(x) of (fH|) can be at most bounded on K. This leads us to 
the following definition. 

Definition 3. We say fi € spec (d x C[u\) is there exists a non-trivial bounded function ip such that 
d x /L[u]ip — flip or, equivalently, if there exists a A € C with |A| = 1 such that 

det ( M(/i) - AI) = 0. 

Following Gardner ^TT$ , we define the periodic Evans function D : C x C — * C to be 

D(fi, A) = det ( M{p) - AI) . 

Moreover, we say the periodic solution u(x;a,E,c) is spectrally stable if spec (d x C[u\) is a subset of the 
imaginary axis. 

Remark 4. In general, one defines a solution u(x;a, E,c) to be spectrally stable if spec(d x £[u]) does not 
intersect the open right half plane. However, the Hamiltonian nature of the problem implies the spectrum 
is symmetric with respect to the real and imaginary axis, from which our definition follows. Notice that 
since we primarily concerned with roots of D(fj,,X) with A on the unit circle, we will frequently work with 
the function D(/j,, e lK ) for (tel. However, the more general definition given above is useful for certain 
analyticity arguments. 

It follows that the set spec(d x £[u]) consists of precisely the L°° eigenvalues of the linearized operator 
d x C[u]. Moreover, if we define a projection operator tk\ : Cx C — > C by 7Ti(zi, Z2) = z\ for {z\,z-z) € C x C, 
then the projection of the zero set of D{[i,e %K ) in CxS 1 via 7Ti is precisely the set spec(i9 a ;£[u]). In 
particular we see that D(/i,l) detects spectra which corresponds to perturbations which are periodic with 
the same period as the underlying solution u. If the zero set of -D(/i, 1) is a subset of the imaginary axis, 
it follows that the underlying periodic solution is spectrally stable with respect to perturbations of the 
same period. Although this is a long way from concluding any sense of nonlinear stability, we find this 
analysis is vital to understanding the orbital stability calculation. In particular, we recall the following 
easily proved result. 

Lemma 1. The function D{fi, 1) is odd and the limit lim^^oo sign (-D(/i, 1)) exists and is negative. 

The proof relies on an asymptotic formula for the quantity tr( M(/i)) for /i» 1, as well as a structural 
property of the Evans function which holds due to the Hamiltonian nature of the problem. For details, 
see Proposition 1 and Lemma 1 of [S]. In order to determine if D(/j,, 1) has non-zero real roots, we also 
need the following technical lemma. 

Lemma 2. The periodic Evans function satisfies the following asymptotic relation near /i = 0: 

D(ji, 1) = -\{T, M, P} Q , £ , C M 3 + 0(\fi\ 4 ). 
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The proof of this lemma is somewhat technical, although straightforward, and will not be repro- 
duced here: see Theorem 3 in [8] for details. With the above two lemmas, it immediately follows that 
if {T, M, P} a ,E,c is negative, then the number of positive roots of -D(/x, 1) is odd and hence one has 
exponential instability of the underlying periodic traveling wave. Moreover, we will show in Lemma 2] 
that Te > implies C[u] has exactly one negative eigenvalue. It follows from that the linearized operator 
d x C[u] has at most one unstable eigenvalue with positive real part, counting multiplicities (see Theorem 
3.1 of [20j). Since the spectrum of d x C[u] is symmetric about the real and imaginary axis, it follows 
that all unstable periodic eigenvalues of the linearized operator must be real. This proves the following 
extension of Corollary 1 in [S] . 

Theorem 2. Let u(x;ao, Eq,co) be aperiodic traveling wave solution of ([T]). If Te(o,q, Eq, Co) > 0, then 
the solution is spectrally stable to perturbations of the same period if and only if {T, M, P} a ,E,c is positive 
at (a ,E Q ,c ). 

Notice that if Tg < the operator C[u] has two negative eigenvalues and hence if {T, M, P} a ,E,c > 
there is no way of proving from these methods whether the number of periodic eigenvalues of d x C\u] with 
positive real part is equal to zero or two. In the next section, we will prove that, if {T, M, P} a ,E,c is 
positive, then the periodic traveling wave u(x; a, E, c) is orbitally stable if and only if the sub-determinants 
Te and {T, M} aj £ are also positive. We will show these sub-determinants are positive in certain cases, 
but we do not know if this is true in general. We conjecture, however, that this condition holds for all 
(a, E, c) € fi in the case of a power-law nonlinearity. 



4 Orbital Stability 

In this section, we prove our main theorem on the orbital stability of periodic traveling wave solutions of 
(fTj). Throughout this section, we assume we have a T-periodic traveling wave solution u(x; do, Eq, cq) of 
equation (JTJ) , i.e. we assume u satisfies 

l -u 2 x + F(u) - °^u 2 - a u = E (15) 

with (cio, -Eci; Co) € ^ and T = T(ao, Eq, cq). Moreover, we assume the non-linearity / present in ([T]) is 
such that the Cauchy problem for is globally well-posed in a real Hilbert space X of real valued T 
periodic functions defined on K, which we equip with the standard L 2 ([0,T]) inner product 

id, h) := / g(x)h(x)dx 
Jo 

for all g,h £ X, and we identify the dual space X* through the pairing 

= / g{x)h(x)dx 
Jo 

for all g € X* and h € X. In particular, notice that L 2 ([0,T]) is required to be a subspace of X. For 
example, if f(u) = u 3 /3, corresponding to the modified Korteweg-de Vries equation, then the Cauchy 
problem for Q is globally well-posed in the space 

ff p s er ([0,T];E) = { 3 6F([0,T];1) : g(x + T) = g(x) a.e.} 

for all s > 5, where we identify the dual space with /Z^" e *([0, T];R) through the above pairing (see [5] for 
proof). Moreover, due to the structure of the gKdV, we make the natural assumption that the evolution 
of ([¥]) in the space X is invariant under a one parameter group of isometries G corresponding to spatial 
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translation. Thus, G can be identified with the real line K acting on the space X through the unitary 
representation 

(R e g)(x) = g(x + $ 

for all g £ X and £ £ G. Since the details of our proof works regardless of the form of the non-linearity 
/, we make the above additional assumptions on the nonlinearity and make no other references to the 
exact structure of the space X nor /. 

In view of the symmetry group G, we now describe precisely what we mean by orbital stability. We 
define the G-orbit generated by u to be 

O u := {R^u : £ £ G}. 

Now, suppose we have initial data cf>o £ X which is close to the orbit O e . By orbital stability, we mean 
that if <p(-,t) £ X is the unique solution with initial data 4>o, then </>(•, t) is close to the orbit of u for all 
t > 0. More precisely, we introduce a semi-distance p defined on the space X by 

p{g,h) = mi \\g- Rf.h\\ x , 

and use this to define an e-neighborhood of the orbit O u by 

U e := {4> £ X : p(u, 4>) < s}. 
The main result of this section is the following restatement of Theorem ([I]) . 

Proposition 1. Let u(x) — u(x;ao,Eo,co) solve (|15l) and suppose the quantities Te, {T, M} a E, and 
{T,M,P} a E,c o,re positive at (ao, Eq,cq). Then there exists positive constants Co,£o such that if <f>o £ X 
satisfies p(4>o,u) < e for some e < Eq, then the solution <fi(x,t) of ([T]) with initial data 4>o satisfies 
p(<f>(;t),u) < C e. 

Remark 5. Notice that Theorem^ implies a periodic solution u(x; ao, Eq, Cq) of ([5]) is an exponentially 
unstable solution of ([l]) if {T, M, P} a ,E,c is negative at (ao, E$, cq). Thus, the positivity of this Jacobian 
is a necessary condition for nonlinear stability. 

Before we prove Proposition Q] we wish to shed some light on the hypotheses. Recall that the classical 
action K(a, E, c) of the periodic traveling wave satisfies 

V a ,E, c K{a, E, c) = \ M(a, E, c),T{a, E, c), ^P(a, E, c) 

As a result, we can write its Hessian as 

Dl^ c K{a,E,c) = 

Proposition [1] thus states that if (ao, Eq,cq) £ £1, the corresponding periodic traveling wave solutions 
of ([T]) is orbitally stable if the principle minor determinants of D\ E c K(a, E, c) satisfy di = Te > 0, 
c?2 = {M, T} a ^E < 0, and c?3 = {M, T, P} a ,E,c < 0. It is clear that a necessary condition for this claim 
is that the Hessian D\ E c K(a, E, c) is invertible with precisely one negative eigenvalue. However, this is 
clearly not sufficient. Although we are unable to give a complete characterization of the orbital stability 
in terms of orthogonal invariants of D\ E c K(a, E, c), such a connection would indeed be very interesting. 

We now proceed with the proof of Proposition[TJ which follows the general method of Bona, Souganidis 
and Strauss [6], and consequently that of Grillakis, Shatah and Strauss |13 J .|15 j . We define the following 
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functionals on the space X, which correspond to the "energy", "mass" and "momentum" respectively: 

£{<t>) ■= £ Q^(«) a -*W*))) dx 

M{<p) := [ 4>(x) dx 
Jo 

These functionals represent conserved quantities of the flow generated by |T]). In particular, if <p(x, t) is a 
solution of H]) of period T, then the quantities £((f>(-, t)), M. (</>(•,£)), and V(4>{-,t)) are constants in time. 
Moreover, notice that £(u) = H(ao,E ,co), M(u) — M (a , E , cq), and V(u) = P(a , E ,c ) where H, 
M and P are defined in (P|)- (fTTj) . 

Remark 6. Throughout the remainder of this paper, the symbols M and P will denote the functionals 
M. and V restricted to the manifold of periodic traveling wave solutions of |T]) with (a, E, c) G fi. 

Remark 7. Calculations in similar vain have been carried out recently in the special cases of cnoidal 
solutions of the KdV JBjj, as well as for traveling wave solutions of the modified KdV arising from |T]) 
with f(u) = u 3 |7J/. In each of these cases, however, it was assumed that a = 0, or equivalently that 
M(a,E,c) = 0. While this is always possible for the KdV (due to Galilean invariance), this is not 
possible for general nonlinearities without restricting your admissible class of traveling wave solutions, 
i.e. restricting Q. As we are interested in deriving universal conditions for stability of traveling wave 
solutions of ([T]), we are forced to work with all three functionals defined above. 

It is easily verified that £, M. and V are smooth functionals on X, whose first derivatives are smooth 
maps from X to X* defined by 

£ \<f) = - /(</>), M'(4>) = 1, V'{4>) = 4>- 

If we now define an augmented energy functional on the space X by 

£ (<t>) := £{cfy + c V((b) + a M{4>) + E T (16) 
it follows from (fT5f that £q(u) — and £' (u) — 0. Hence, u is a critical point of the functional £q. 

Remark 8. Notice that the added factor of EqT on the right hand side of (|16[) is not technically needed 
for our calculation. However, we point out that (formally) if we consider variations in £q in the period 
we obtain 

^oW>)|, =u = \u 2 AT) - F(u(T)) + au(T) + E+ (^(u), ^ 

= \ul{T) + \ul(T) + (£'{u),^ 
= 

since u x (T) = and £q(u) — 0. Hence u is also (formally) a critical point of the modified energy with 
respect to variations in the period. It would be very interesting to try to make this calculation rigorous 
and to see if it allows one to extend orbital stability results to include perturbations with period close to 
the period of the underlying periodic wave. Again, this is clearly all very formal and we will make no 
attempt at such a theory here. 
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To determine the nature of this critical point, we consider its second derivative £(f , which is a smooth 
map from X to C(X,X*) defined by 

£q(4>) = -<t>xx - /'(</>) + c - 

This formula immediately follows by noticing the second derivatives of the mass, momentum, and energy 
functionals are smooth maps from X to C(X,X*) defined by 

£"(</>) =-%-f'{<f>), M"(0) = O, P"(0) = 1. 

In particular, notice the second derivative of the augmented energy functional £q at the critical point u 
is precisely linear operator C[u] arising from linearizing © with wave speed cq about u. It follows from 
the comments in the previous section that £q(u) is a self-adjoint linear operator on Lp Cr ([0, T]; K) with 
compact resolvent. In order to classify u as a critical point of £q, we must understand the nature of the 
spectrum of the second variation C[u\: in particular, we need to know the number of negative eigenvalues. 
This is handled in the following lemma. 

Lemma 3. The spectrum of the operator C[u] considered on the space Lp er ([0,T]) satisfies the following 
trichotomy: 

(i) If Te > 0, then C[u] has exactly one negative eigenvalue, a simple eigenvalue at zero, and the rest 
of the spectrum is strictly positive and bounded away from zero. 

(ii) If Te — 0, then C[u] has exactly one negative eigenvalue, a double eigenvalue at zero, and the rest 
of the spectrum is strictly positive and bounded away from zero. 

(Hi) IfTs < 0, then C[u] has exactly two negative eigenvalues, a simple eigenvalue at zero, and the rest 
of the spectrum is strictly positive and bounded away from zero. 

Proof. This is essentially a consequence of the translation invariance of |T]) and the Strum-Liouville 
oscillation theorem. Indeed, notice that for any (eG the function R^u is a stationary solution of (JBJ) 
with wave speed Co and a — clq. Differentiating this relation with respect to £ and evaluating at £ = 
implies that £[ii]u :I = 0. Moreover, since u is radially increasing on [0, T] from its local minimum there, 
u x is periodic with the same period as u and hence u x S Lp Cr ([0, T]). This proves that zero is always a 
periodic eigenvalue of C[u] as claimed. To see there is exactly one negative eigenvalue, notice that since 
u is T-periodic with precisely one local critical point on (0,T), its derivative u x must have precisely one 
sign change over its period. By standard Strum-Liouville theory applied to the periodic problem (see 
Theorem 2.14 in [T7]), it follows that zero must be the either the second or thircH eigenvalue of C[u] 
considered on the space Lp er (R). 

Next, we show that zero is a simple eigenvalue of C[u] on the space Lp Or ([0, T]) if and only if Te ^ 0. 
To this end, notice that the periodic traveling wave solutions of (© are invariant under changes in the 
energy parameter E associated with the Hamiltonian ODE ©. As above, it follows that C[u]ue = 0. 
We must now determine whether the function ue belongs to the space Lp Or ([0, T]). Since it is clearly 
smooth, we must only check whether it is periodic with the same period as the underlying wave u. To 
this end, we use u x and ue as a basis to compute the monodromy matrix m(0) corresponding to the 
equation C[u]v = 0. Notice that differentiating the relation E = V(u-;a,c) with respect to E and 
evaluating at (ao, E$, cq) gives ^^-V'(u^; ao, Co) = 1, and hence ^jjr is non-zero at (ao, Eq, cq). Defining 

yi{x) = (^jj^j ue and yzix) — — (V'(u— ; do, Co)) 1 u x (x), it follows from direct calculation that 

yx(0) = l, 2/ 2 (0) = 0, 
yj(0) = 0, 2/2(0) = 1. 



5 Clearly, we mean with respect to the natural ordering on R. 
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Thus, it follows by calculating ue(T) by the chain rule that we have 

■»<»> = ( I ?)• 

where again we have used the fact that V'(u-; a, c)^p§- = 1- Thus, it follows that zero is a simple 
eigenvalue of C[u] if and only if Te ^ 0, and that the multiplicity will be two in the case Te = 0. 

Finally, to determine whether fj, = is the second or third eigenvalue of C[u], we note that from the 
results of [5] we have 

sign(T B )= S igii(tr(m / ,(0))). (17) 

By the oscillation theorem for Hill-operators with periodic coefficients (see Theorem 2.1 in ^TT\), it follows 
that \i — is the second periodic eigenvalue of C[u] if and only if Te > 0, and the third periodic eigenvalue 
if and only if Te < 0. This completes the proof. □ 

Notice that in the solitary wave case, the spectrum of the operator C[u] always satisfies (i) in the 
above trichotomy. Since E is not restricted to be zero in the periodic context, it is not surprising that 
such a non-trivial trichotomy might exist. The next lemma shows that for a large class of nonlinearities, 
the period is indeed an increasing function of E within the region f2. 

Lemma 4. Let (ao, Eq,cq) £ O and u — u( ■ ; ao, Eq, cq) denote the corresponding periodic solution of 
(|5|) with wave speed cq and period T = T(ao, Eq, cq). If the nonlinearity f in ([![) is such that f'(u) is 
co-periodic with u, then Te > at (ao,i?Oj c o)- 

Proof. If f'(u) is co-periodic with u, then the operator C[u] is a Hill operator with even potential with 
period T = T(ao, Eq, cq). Thus, u x is a periodic eigenvalue of C[u] which satisfies Dirichlet boundary 
conditions. Moreover, since u x changes signs once over (0, T), it follows that zero is either the second or 
third periodic eigenvalue of C[u\. Since the first periodic eigenvalue must be even, and hence must satisfy 
Neumann boundary conditions, it follows by Dirichlct-Neumann bracketing |21j that zero must be the 
second periodic eigenvalue of C[u]. Using the notation of Lemma [3J it follows that tr( m M (0)) > and 
hence T E > by equation |T7]). □ 

Remark 9. In particular, it follows that if f(u) = u p+1 for some p > 1 then the spectrum of C[u] will 
satisfy (i) in Lemma\^ 

Throughout the rest of the paper, unless otherwise stated, we will assume that Te > at (ao, Eq, cq) 
and hence zero is a simple eigenvalue of the operator C[u] considered on the space Lp Cr (R). In particular, 
we assume that the map E — > T(E, ao,Ca) does not have a critical point at Eq. It follows that we can 
define the spectral projections n_, Ho and n + onto the negative, zero, and positive subspaces of the 
operator C[u] (respectively) via the Dunford calculus. Thus, any <p S X can be decomposed as a linear 
combination of u x , an element in the positive subspace of C[u], and x, where x is the unique positive 
eigenfunction of C[u] with ||x||l 2 ([o.t]) = 1 which satisfies 

(£MX,X) = -A 2 

for some A > 0. From the above definition of \ it follows that x is the eigenfunction corresponding to 
the unique negative eigenvalue —A 2 of C[u]. 

From Lemma [H we know that u is a degenerate saddle point of the functional £q on X, with one 
unstable direction and one neutral direction. In order to get rid of the unstable direction, we note that 
the evolution of ([1]) does not occur on the entire space X, but on the co-dimension two submanifold 
defined by 

E := G X : M{<j>) = M(a , Eq, cq), V(4>) = P(a , Eq, c )}. 
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It is clear that So is indeed a smooth submanifold of X in a neighborhood of the group orbit O u . 
Moreover, the entire orbit O u is contained in So- The main technical result needed for this section is 
that the functional £$ is coercive on So with respect to the semi-distance p, which is the content of the 
following proposition. 

Proposition 2. If each of the quantities Te, {T, M} a ^E, and {T, M, P} a ,E,c are positive, then there 
exists positive constants C\, S which depend on (ao, Eq,cq) such that 

SoW-W^dp^u) 2 

for all <f> € So such that p(<j),u) < 5. 

The proof of Proposition's broken down into three lemmas which analyze the quadratic form induced 
by the self adjoint operator C[u]. To begin, we define a function <pQ by 

<f> (x) := {u(x\a,E,c),T(a,E,c),M(a,E,c)}a iE ,c\( a0tEo>Co y 

It follows from a straightforward calculation (see Proposition 4 in [5]) that <f>o £ X and 

£k# = {T, M} E ,c - {T, M} a , B u, 

where the right hand side is evaluated at (ao, Eq, en). This function plays a large role in the spectral 
stability theory for periodic traveling wave solutionqj of (fTJ) outlined in section 3. In particular, we have 
d x £[u]<fio — —{T, M} a ,EU x , and hence, assuming {T, M} a> E ^ at (ao, E$, cq), 4>q is in the generalized 
periodic null space of the linearized operator d x C[u]. See [H] for more details on the role of this function 
in the spectral stability theory outlined in section 2. 
Notice that <fio does not belong to the set 

% = {<t>eX:(u^) = (1,^ = 0}, 

which is precisely the tangent space in X to So. Indeed, while (l,4>o) = by construction, the inner 
product (u, 4>q) — {T, M, P} a ,E.c does not vanish by hypothesis. Using the spectral resolution of the 
operator C[u], we begin the proof of Proposition [5] with the following lemma. 

Lemma 5. Assume that Te > and {T, M} Ut E > 0. Then 

(£M<M) >0. 

for every <f> S Tq which is orthogonal to the periodic null space of C[u\. 

Proof. The proof is essentially found in [6]. First, suppose that Te > and note that by Lemma [3] we 
can write 

(p = a\ + [3u x + p 
(j) = Ax+p 

for some constants a, j3, A, and functions p and p belonging to the positive subspace of C[u]. By 
assumption the quantity 

(£k#o, 0o) = ~{T, M} a , E {T, M, P} a ^ c (18) 
is negative, and hence the above decomposition of 4>o implies that 

> (-\ 2 ax + C[u]p,ax + l3u x +p) = —A 2 a 2 + (C[u]p,p) , (19) 

^Actually, the function u c plays a large role in our analysis via the periodic Evans function. However, since u c is not in 
general T-periodic due to the dependence of the period on the wave speed, we work here with its periodic analogue rf>o . 
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which gives an upper bound on the positive number (C[u]p,p). Similarly, using the above decomposition 
of <fi we have 

= {£[u](po, 0) = -X 2 Aa + {£[u]p,p) . (20) 
Therefore, a simple application of Cauchy-Schwarz implies 

(£N<M =-X 2 A 2 + (C[u]p 7 p) 

„ „ (X 2 aA) 2 
> —X A + V x9 / 

= 

as claimed. 

In the case where Te — 0, there exists a ip G X linearly independent from u x such that C[u]i]j = 0, 
and thus </>o an d 4> admit representations 

4>a = ax + + P u x + P 
<j) = Ax + p 

for some constant 71. It is easily verified that equations (fTO)) and (|20p still hold, and hence 

(£[<#, 0) = (-X 2 A X + C[u]p, A X +P) > 
as before. □ 

Our strategy in proving Proposition [3 is to find a particular set of translates of a given (j> E U £ for 
which the inequality holds. To this end, we find a set of translates for each <j) £ U e such that this set is 
orthogonal to the periodic null space of the second variation C[u] of the augmented energy functional. 
This is the content of the following lemma. 

Lemma 6. There exists an e > and a unique C 1 map a : U £ — > K such that for all <fi € U £ , the function 
</>(■ + o:((f))) is orthogonal to u x . 

The proof is presented in [6], and is an easy result of the implicit function theorem. Indeed, if we 
define the functional 77 : X x K — > R by 

f T 

T)((f),a)= / 4>(x + a)u x {x)dx, 
Jo 

then ^77(0, Qf)|/^ )— ( u 0) = u ^ a ' ^ and hence the lemma follows by the implicit function theorem 
and the fact that by translation invariance the function a can be uniquely extended to lA e for e > 
sufficiently small. 

We now complete the proof of Proposition [2] by proving the following lemma. 

Lemma 7. If each of the quantities Te, {T, M} a ^E, and {T, M, P} a ,E,c o,re positive, there exists positive 
constants C and e such that 

£ { < j>)-£ {u)>CU{- + a {4>))-u\\ 2 x 

for all 4>eU e n E Q . 

Proof. Let e > be small enough such that Lemma O holds. Fix </> e W E fl Eo and write 
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where y G Tq. Moreover, define v — cf>(- + «(</>)) — u and note that by replacing u with R^u if necessary 
we can assume that \\v\\x < s. By Taylors theorem, we have 

M (oo, E , c ) = M(<f>) = M(a , So, c ) + (1, v) + O {pfx) ■ 

Since (1, v) = (3T, it follows that /3 = O (\\v\\ x ). Similarly, we have 

P(ao, E , c ) = P(a , E ,c ) + {u, v) + O (\\vf x ) . 

Moreover, a direct calculation yields 



(u,v) = 7 ||u|||2 ([0]T]) - 



T 



■p(u), 



where (u) = u(x)dx. Since (u) 2 < T\\u\\^ 2 ^ T ^ by Jensen's inequality, it follows that 7 = (IHIx)- 
Now, by Taylor's theorem and the translation invariance of £q, we have 

foW = So (<M- + a{4>))) 

= £ (u) + (£q («),«) + i (^>K u) + o (\\v\\ 2 x ) 
= £ (u) + ~(£[u}v,v)+o(\\v\\ x ). 
Hence, by the previous estimates on 7 and (3, it follows that 

Eotf) - £o(u) = \ (C[u]v, v) + o (\\v\\ 2 x ) 
= ±(£y,y)+o(\\v\\ x ). 
Since y £ Tq and {y,u x ) = by Lemma H3 it follows from Lemma [5] that 

£o(J>)~£o(u)> ^||y|| 2 + o (\\v\\ 2 x ) . 



Finally, the estimates 



\y\\x 



v — ju — (3 - 



T 



> 



\v\\x 



71* — - 



X 

7(u) 



T 



x 



> \\v\\ x -o 



X 



prove that £q(4>) — £q(u) > ^\\v\\ x for \\v\\x sufficiently small, and hence completes the proof. 



□ 



Proposition [2] now clearly follows by Lemma [7| and the definition of the semi-distance p. It is now 
straightforward to complete the proof of Proposition [TJ 

Proof of Proposition We now deviate from the methods of [5] , [T3] and [T5] , and rather follow the 
direct method of [10]. Let 5 > be such that Proposition [2] holds, and let e € (0, S). Assume c/)q G X 
satisfies p(<j>o,u) < e for some small e > 0. By replacing </>o with R^4>q if needed, we may assume that 
||0O — u \\x < £■ Since u is a critical point of the functional £q, it is clear that we have 



£b(to) - £b(«) < C x e 2 
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for some positive constant G\. Now, notice that if 4>o € ^o> then the unique solution 4>{-,t) of (fT]) with 
initial data <fio satisfies </>(■, t) £ So f° r ah i > 0. Thus, Proposition [2] implies there exists a C2 > such 
that p((f>(- , t) , u) < C2S for all i > 0. Thus </>(■, i) € W e for all i > 0, which proves Proposition Q] in this 
case. 

If 4>o £ Xjq , then we claim we can vary the constants (a, E, c) slightly in order to effectively reduce this 
case to the previous one. Indeed, notice that since we have assumed {T, M, P} a ,E,c 7^ at (do, Eq, cq), it 
follows that the map 

(a, E, c) 1 ^ (T(u( ■ ; a, £, c)), M (u( ■ ; a, £?, c)), P(u( ■ ; a, E, c))) 

is a diffeomorphism from a neighborhood of (do, £t)) Co) onto a neighborhood of (T, M(ao, -Eoj Co), -P(dO) ^O) Co 
In particular, we can find constants a, -E, and c with |a| + |£"| + \c\ = 0(e) such that the function 

u = u( ■ ; a + a, Eq + E , c + c) 

solves ([T]), belongs to the space X, and satisfies 

M(a + a, E + E, c + a) = M{<j) ) 
P{a Q + a,E + E, c + c) = V{(j) ). 

Defining a new augmented energy functional on X by 

£{<t>) = £ (4>) + cV(4>) + aM{(f>) + ET, 

it follows as before that 

for some C3 > as long as p{4>{-, t),u) is is sufficiently small. Since u is a critical point of the functional 
£ we have 

C 3P (<X-,i),u) 2 < £(<h) " £ (S) < c 4 Ho - ufx 
for some C4 > 0. Moreover, it follows by the triangle inequality that 

||0o - u\\x < ||0o - «||x + II" - u\\ x < C 5 e 
for some C5 > and hence there is a C§ > such that 

p{<j){-,t),u) < p((p(-,t),u) + \\u-u\\ x < C* 6 £ 

for all t > 0. The proof of Proposition [1] and hence Theorem Q] is now complete. □ 

We would like to point out an interesting artifact of the above proof. Notice that the step at which 
the sign of the quantities {T, M, P} a ,E,c and {T,M} a> E came into play was in the proof of Lemma O 
from which we have the following corollary. 

Corollary 1. On the set Q, the quantity {T, M} a ,£ is positive whenever {T, M, P} a ,E,c is negative and 
Te is positive. 

Proof. This is an easy consequence of Theorem[5]and equation (TT5|) . Indeed, if {T, M} a> E and {T, M, P} a ,E,c 
were both negative for some (ao,-Et)7 c o) G n, then by the proof of Lemma [5] we could conclude that 
• ; do, Eq, cq)]4>, 4>) > for all (f> € 7o which are orthogonal to u x . Since this is the only time in 
which the signs of these quantities arise, it follows that Proposition [T] would hold thus contradicting 
Theorem ©. □ 



5 ANALYSIS NEAR HOMOCLINIC AND EQ UILIBRIUM SOL UTIONS 



17 



It follows that we have a geometric theory of the orbital stability of periodic traveling wave solutions 
of |lj to perturbations of the same period as the underlying periodic wave. In the next two sections, 
we consider specific examples and limiting cases where the signs of these quantities can be calculated. 
First, we consider periodic traveling wave solutions sufficiently close to an equilibrium solution (a local 
minimum of the effective potential) or to the bounding homoclinic orbit (the separatrix solution). By 
considering power-law nonlinearities in each of these cases, we give necessary and sufficient^ conditions for 
the orbital stability of such solutions. Secondly, we consider the KdV and prove that all periodic traveling 
wave solutions are orbitally stable to perturbations of the same periodic as the underlying periodic wave. 

5 Analysis Near Homoclinic and Equilibrium Solutions 

In this section, we use the theory from section 4 in order to prove general results about the stability of 
periodic traveling wave solutions of ([T]) in two distinguished limits: as one approaches the solitary wave, 
i.e. (a, E,c) £ Cl and consider the limit T(a, E, c) — > oo for fixed a, c, as well as in a neighborhood of 
the equilibrium solution, i.e. near a non-degenerate local minimum of the effective potential V(u;a,c). 
Throughout this section, we will consider only power-law nonlinearities. 

We begin with considering stability near the solitary wave. Our main result in this limit is that the 
quantities Te and {T, M} Qj £ are positive for (ag, Eq^ cq) £ Cl with sufficiently large period. Hence, the 
orbital stability of such a solution in this limit is determined completely by the periodic spectral stability 
index {T, M, P} a .E.c, which in turn is controlled by the sign of the solitary wave stability index This 
is the content of the following theorem. 

Theorem 3. In the case of a power-law nonlinearity, i.e. f(u) = u p+1 with p > 1, a periodic trav- 
eling wave solution of §Q of sufficiently large period and (a,E,c) £ Cl is orbitally stable if p < 4 and 
exponentially unstable to perturbations of the same period as the underlying wave if p > 4. 

Proof. By the work in [8], the quantity M a is negative for such (a,E,c) £ Cl. Moreover, since we are 
working with a power-law nonlinearity, the periodic traveling wave solutions satisfy the scaling relation 



as CI 9 (a, E, c) — > (0, 0, c) for a fixed wave speed. Since {T, M} a> E = M E — TsMa, it follows from Lemma 
|4]and Theorem [2] that the solutions u(x; a, E, c) with (a, E,c) £ Cl of sufficiently large period are orbitally 
stable if p < 4 and exponentially unstable to periodic perturbations if p > 4. See [5] for more details of 



Next, we consider periodic traveling wave solutions near the equilibrium solution. We will use the 
methods of this paper to prove that such solutions are orbitally stable to periodic perturbations, provided 
that a is sufficiently small. To begin, we fix a wave speed c > and consider (|TJ) with a power-law 
nonlinearity f(u) — u p+1 with p > 1. Recall that Te > by Lemma IH and hence it suffices to prove 
that {T, M} a ^E and {T, M, P} a ,E,c are both positive near the equilibrium solution. By continuity, it is 
enough to evaluate both these indices at the equilibrium and to show they are both positive there. This 
is the content of the following lemma. 

7 Exccpt in the exceptional case of being near the homoclinic orbit for p = 4. 




from which we get the asymptotic relation 




these asymptotic calculations. 



□ 
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Lemma 8. Consider ([T]) with a power-law nonlinearity f{u) = u p+l for p > 1. Then the quantity 
M a is negative for all (oq, Eq, Co) G f2 sucft i/iai \a\ is sufficiently small and the corresponding solution 
u{ ■ ;ao, Eq,cq) is sufficiently close to the equilibrium solutio^. 

Proof. First, denote the equilibrium solution as u a c and let E*(a, c) = V{u a>c ; a, c). It follows that 

lim T(a, E, c) = —— 
e\e* ^/bp 

and that the equilibrium solution admits the expansion 



1 + - +0{a 2 ). 
Pj 

Now, solutions near the equilibrium u a ^ c can each be written as u(x;a, E,c) = P a ,E,c{k a ,E,cx), where 
k a ,E,cT(a, E, c) = 2n and P a ,E,c is a 2tt periodic solution of the ordinary differential equation 

k 2 a ,E,cv" + v p+1 - c 1+1/p a = 

such that 

Pa,E*,c = Uo,c, k 2 a E ^ c = {p+ l)u P a c - 1. 

Straightforward computations give the expansions 



P a ,£,cO) = u a , c + 0{^E-E* (1 + a 2 )), 

k 2 aEc = cp+(p+ l)ca + 0((E - E*) + a 2 ). 



Thus, the mass M(a, E, c) can be expanded as 



r2ir/k a , E ,a 

M(a,E,c)= / Pa,EA k a,E,cZ)dz 



i r 2 * ~ 

, / P a ,E,c( Z ) dz 

Ka,E,c JO 



= ^(l + ^^)+0(^E-EF + a 2 ) 



It follows that 

#M(a,E lC )| mp , 
9a v 'l(o,B*,c) 

which is negative for p > 1 . 

The case p = 1, which corresponds to the KdV equation, will be discussed in the next section. There 
we will show that although M a vanishes at the equilibrium solution, it is indeed negative for nearby 
periodic traveling waves with the same wave speed c, i.e. 

d 2 

-M(a,E,c)\ (n R , . <0. 



dEda K ' ' ;| (o,£*,c) 

□ 

Next, we must determine the sign of the periodic spectral stability index {T, M, P} a ,E,c- Although 
it follows from Theorem 4.4 in [H] that this index must be positive we present an independent proof 

8 In the case of the KdV (p = 1), M a is negative in a deleted neighborhood of the equilibrium solution. 
9 They prove the spectrum of the linearized operator d x C[u] intersects the real axis only at fj, = for such small amplitude 
solutions. 
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based on the periodic Evans function methods of [8 . To this end, we point out that the Hamiltonian 
structure of the linearized operator 9 x £[n] we have the identity 



{T,M,P} a . E . c = --tv(M^(0)) 



where M(/x) is the corresponding monodromy operator (see Theorem 3 of [5] for details). Thus, it is 
sufficient to show that tr( M MAlAl (0)) is negative near the equilibrium solution. This is the content of the 
next lemma. 

Lemma 9. Consider |T]) with a power-law nonlinearity f(u) = u p+1 with p > 1. Suppose uo is a non- 
degenerate local minima of the corresponding effective potential V(u;a,c). Then tr ( M MW (0)) < at 

Uq. 

Proof. The key point is that the matrix H(x, fi) in (|14[) reduces to the constant coefficient matrix 

/ 1 

H(/i) =10 1 

V -H -V"(u ;a,c) 

at the equilibrium solution uq. Thus, the corresponding monodromy operator at uq can be expressed as 
M(/i) = exp (H(/x)To), where To = ^jg- Thus, in order to calculate the function tr(M(/i)), it is sufficient 
to calculate the eigenvalues of the constant matrix H(/z). 

Now, the periodic Evans function corresponding to the constant coefficient system induced by H(/i) 
can be written as 

D (fi, A) = dct ( H(/i) - AI) = -A 3 - V"(u ; a, c)A - \i. 

In particular, notice that J^Do(n, A) = —A 2 — V"(uq\ a, c). Since V"(uo;a, c) > it follows that the 
function Do(fi. ■ ) will have precisely one real root for each /igl. This distinguished root is given by 
the formula 



m 1/3 F"M ( (^ + V3^/27v 2 + W"(u )A 

7i(M)- " 1 



1/3 > 



(9/i + W 2 V + 4^"(uo) : 



1/3 2V332/3 



=:a(jj,) =:/3(m) 

Defining lu = exp(2-7ri/3) to be the principle third root of unity, the two complex eigenvalues of H(/i) 
can be written as 72 (m) = ua(ji) + uJ/3(/x) and 73 (pi) = u7a(/x) + o;/3(/x), and hence 

tr (Mfci)) = exp ( 7 i(/i)T ) + exp (72(^0) + exp (73(^0) . 

Now, a straightforward, yet tedious, calculation using the facts that 1 + uj + 57 = and cj 2 = ZU implies 
that 

tr ( M w/i (0)) = 9T 2 (a"(0)/3'(0) + a/(0)/3"(0)) + 3T 3 (a'(0) 3 + /3'(0) 3 ) . 
Moreover, from the definitions of a and (3 we have 

«' <0 > " -wk^) " ,3 ' (0) ' a " d a " (0) " iF^FTi ^ -""(")• 

Therefore, we have the equality 

tr((M mi (0)) = - (, ' V " 



y"M 9 / 2 ' 

which is clearly negative. □ 
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Remark 10. Notice that the proof of Lemma\$\ holds verbatim if one considers more general nonlinear- 
ities. In particular, one could obtain a similar expression for tr( M w (0)) and show the modulational 
instability index from always vanishes at the equilibrium solutions of the traveling wave ordinary 
differential equation^ Moreover, this proof holds regardless of the size of the parameter a. 

Therefore, it follows that in the case of a power-nonlinearity and solutions sufficiently close to a non- 
degenerate minima of the effective potential, each of the quantities Te, {T, M} a ^E, and {T, M, P} a ,B,c 
are all positive. Therefore, Theorem [T] immediately yields the following result. 

Theorem 4. Consider equation JT]) with a power-law nonlinearity f(u) — u p+1 for p > 1. Then the 
periodic traveling wave solutions u(x; a, E, c) with (a, E, c) £ Q, and a 2 + (E — E*) 2 sufficiently small are 
orbitally stable in the sense of Theorem^ 



6 The Korteweg-de Vries Equation 

In this section, we will apply the general theory from section 4 in order to prove that periodic traveling 
wave solutions of (JT|) with f(u) — u 2 and c > are orbitally stable with respect to periodic perturbations 
if and only if they are spectrally stable to such perturbations. To this end, we notice that solutions of 
the KdV equation 

u t = U X xx + {u 2 ) x - cu x (21) 
are invariant under the scaling transformation 

U,(x, t) I— > c u{\fcx, c 3 ^ 2 t), 

and hence, by ([7]), the stationary periodic traveling wave solutions satisfy the identity 

/ \ / i ii a ^ \ 

u(x; a, E, c) = c u c 1 x; — , — , 1 . 

\ c c A J 

Thus, by scaling we may always assume that c = 1 in (|21| . Moreover, we may always assume that a = 
due to the Galilean invariance of the KdV. Therefore, it is sufficient to determine the stability of periodic 
traveling wave solutions of (|2ip of the form u(x; 0, E, 1). In order to do so, we need the following easily 
proved lemma: 

Lemma 10. Let [i be a (Borel) probability measure on some interval I C R, and let /, g : / — > R be 
bounded and measurable functions. Then 

jf f{x)g[x)d l i - H f(x)d(A (J g(x)d^j = \ j (/(*) ~ /GO) ~ 9(y)) d^ y . (22) 

In particular, if both f and g are strictly increasing or strictly decreasing, and if the support of \x is not 
reduced to a single point, then 

J i f(x)g(x)dfi> nf(x)d(j\ n g (x)di^j . 

The proof of this lemma is a trivial result of Fubini's theorem, as one can see by writing the left hand 
side of (|22p as an iterated integral and simplifying the resulting expression. Now, recall from Lemma U 
that Te > for periodic traveling wave solutions of (|21|1 . To conclude orbital stability, we must identify 
the signs of the Jacobians {T, M} aj E and {T, M, P} a ,E,c- The main technical result we need for this 
section is the following lemma, which uses Lemma [TU] to guarantee the sign of the quantity in (jT5J) is 
completely determined by the Jacobian {T, M, P} a ,E,c- 
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Lemma 11. If f(u) = u 2 in ([T]), then {T, M} a E > for all (ao, Eq, cq) £ £1 which do not correspond to 
the unique equilibrium solution. 

Proof. First, notice that since V a ,E, c K{a 1 E, c) = (M,T,P), it follows that {T,M} a . E = Mf - T E M a , 
and hence by Lemma [Hit is enough to prove that M a < 0. Moreover, by our above remarks it is enough 
to consider the case c = 1 and a = 0. It follows for / given as above, we can find functions u±, U2, u 3 
which depend smoothly on (a, E, c) within the domain il such that 

3(E- V(u; 0, 1)) = (tt - ui)(u - u 2 ){u 3 - u). 

Notice that the assumption that we are not at the equilibrium solution implies that the roots Ui are 
distinct, and moreover that V'{ui\ 0, 1) ^ 0. Since E — V{ui\ 0, 1) = on fi, it follows that 

duj _ Uj 
~da~ ~~ V'(ui; 0,1)' 

Since u\ < and u 2 , u 3 > 0, we have 

du\ „ du2 „ , du 3 „ . „. 

< 0, — ^ < 0, and — i > 0. (23) 

Moreover, since ^1+^2+^3 = ^ we have the relation 

^ + ^=-^>0 (24) 
oa aa oa 

on fi. 

Now, by making the change of variables u i— > s((9) = u 2 cos 2 (0)+it2 sin 2 (6*), we have dii = (u — u 2 )(u 3 
and hence we may express the mass of u(x; a, E, c) as 

rMa u du 



M(a,E,c) = \/2 / 



V-B- V(u;o,c) 



2V6T 72 . (25) 



Notice we suppress the dependence of s(9) on the parameters (a,E,c). Defining a(9) = y/ s(9) — u±, a 
straightforward computation using (f2~4"|) shows that the derivative of the integrand in (f2"5|) with respect 
to the parameter a can be expressed as 

d ( s(9) \ _ du 2 / cos 2 (6>) \ du 3 / sin 2 (6>) 
3a ^3(0) -wj ~ da { 2a(9) J + da { 2a(9) 

(du 2 du 3 \ ( s(9) \ mi (du 2 2 , du 3 . 2 

_ 0U2 / cos 2 (6>)-sin 2 (6>) \ / ou 2 <9m 3 \ / sin 2 (6>) \ _ s(9) 

da \ 2a{9) ) + \ da + da ) \2a(9) 3 ) ~ \ da + da ) 2a(9f 
( du 2 / cos 2 (6>) -sin 2 (6>) \ ( du 2 du 3 \ / sin 2 (6>) \\ 
~ Ul V da V M#) 3 ) + \da + da)\ 2a{9f ) ) ' 

With a little more algebra, this may be rewritten as 

d ( s{9) \ du 2 /cos 2 (6>) -sin 2 (6>)\ du 2 /cos 2 (6>) - sin 2 (60 



9a \y s (0) - Ul J da V 2<r(0) / 9a V 2a(0)3 

'0u a aits \ / s(0) cos 2 ((9) - Ul (cos 2 (0) - sin 2 (0)) ' 



da da J \ 2a(9) s 
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Since the functions cos 2 (#) — sin 2 ((9) and <j(9) 1 arc strictly decreasing on the interval (0, 7r/2), it follows 
from Lemma [TOl that 

for any m > 0. Evaluating the above expression at (a, 25, c) = (0, 25, 1) € fi implies that s(0) > for all 
6 € (0,7r/2), and hence ((Ml) and imply that 




d6> < 



y«(0)-«i, 

at (0, E, 1), from which the lemma follows. □ 

Therefore, our main theorem on the stability of periodic traveling wave solutions of the Korteweg-de 
Vries equation follows by Theorem [TJ Theorem [21 and Lemma [TT1 

Theorem 5. Let u(x; ag, Eq, cq) be a periodic traveling wave solution of (|2ip . i.e. let (ao,Eo,co) G fi. 
Then u is orbitally stable in the sense of Theorem [7] if and only if the solution is spectrally stable to 
perturbations of the same period, i.e. if and only if {T, M, P} a ,E,c > at (do, Eq,Cq). 

An interesting corollary of Theorem [5] applies to cnoidal wave solutions of the KdV. It was suggested 
by Benjamin [3] that such solutions should be stable to perturbations of the same period. This conjecture 
has indeed been proved both by using the complete integrability of the KdV [18] [7] and by variational 
methods as in the present paper [3J. In particular, in [7] it was shown that the cnoidal solutions of the 
KdV are spectrally stable to localized perturbations, and are linearly stable to perturbations with the 
same period as the underlying wave. Clearly then such solutions are spectrally stable with respect to 
periodic perturbations. Paired with Theorem this provides another verification Benjamin's conjecture 
in the case where the cnoidal wave has positive wave speed. 

Corollary 2. The cnoidal wave solutions of (QJ with f(u) = u 2 of the form 

u(x, t) = u + 12k 2 n 2 cn 2 (k (x - x - (8fc 2 K 2 - 4k 2 + u ) t)) , 

with k G [0, 1) and k, Xq, and u$ real constants, are orbitally stable in the sense of Theorem^ if the wave 
8k 2 k 2 — 4k 2 + uq is positive. 



7 Concluding Remarks 

In this paper, we extended the recent results of [S] in order to determine sufficient conditions for the orbital 
stability of the four-parameter family of periodic traveling wave solutions of the generalized Korteweg-de 
Vries equation ([T]). By extending the methods of [6] to the periodic case, a new geometric condition 
was derived in terms of the conserved quantities of the gKdV flow restricted to the manifold of periodic 
traveling wave solution, and it was shown how this could be translated to a condition on the Hessian 
of the classical action of the ordinary differential equation governing the periodic traveling waves. As a 
byproduct of this theory, it was shown that such solutions of the KdV are orbitally stable to perturbations 
of the same period as the underlying wave if and only if they are spectrally stable to periodic perturbations 
of the same period. 

There are several points in which this theory is still lacking. First off, it is not clear what happens 
in the case Tg < 0. In the solitary wave theory, the existence of two negative eigenvalues of the second 
variation C[u] indicates instability. Also, if Te > and {T,M} a E < it is not clear whether this 
implies orbital instability, although we conjecture this is indeed the case. We would like to show in the 
case {T, M} a ^{T, M, P} a ,E,c < and Te > that there exists a 1-parameter family of functions in E 
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which contain the solution u(x; ao, Eg, cq) such that the augmented energy functional £q has a strict local 
maximum at (ao, Eq, cq). However, it is not clear how to do this in a reasonable manner: mainly, one 
must fix the period, mass, and momentum along this curve, and the derivative of this curve at (ao, Eg, cq) 
must also be in X . The existence of such an instability would stand in stark contrast to the solitary wave 
case, where the orbital stability is equivalent to the spectral instability (except possibly on the transition 
curve). However, it seems quite possible that such a situation arises due to the fact that the solitary 
waves are a co-dimension two subset of the family of traveling wave solutions. 

Acknowledgements: The author would like to thank Jared C. Bronski for many useful discussions, 
especially as concerned with the proof of Lemma HI and for looking over the preliminary work regarding 
this paper. 
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